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In open quantum systems, phenomenological master equations with unknown parameters are 
often introduced. Here we propose a time-independent procedure based on quantum tomography 
to reconstruct the potentially unknown parameters of a wide class of Markovian master equations. 
According to our scheme, the system under investigation is initially prepared in a Gaussian state. At 
an arbitrary time t, in order to retrieve the unknown coefficients one needs to measure only a finite 
number (ten at maximum) of points along three time-independent tomograms. Due to the limited 
amount of measurements required, we expect our proposal to be especially suitable for experimental 
implement at ions. 
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I. INTRODUCTION 

Tomographic maps [TJ can be considered a very useful 
tool for reconstructing the physical state or some other 
properties of many physical systems, both in a classical 
(e.g. medical physics, archaeology, biology, geophysics) 
and in a quantum perspective (e.g. photonic states [2], 
photon number distributions [3J SI [S], longitudinal mo- 
tion of neutron wave packets [B]). 

The tomographic analysis is based on a probabilistic 
approach towards physical system investigation. In par- 
ticular, its key ingredient is the Radon transform [TJ. 
Given the phase-space of the system, this invertible inte- 
gral transform allows to retrieve the marginal probability 
densities of the system, i.e. the probability density along 
straight lines. However, while in the classical regime the 
state of the system can be fully described by means of 
a probability distribution on its phase space, this is no 
longer the case of quantum systems. Indeed, due to the 
Heisenberg uncertainty relation, it is not possible to write 
a probability distribution as a function of both momen- 
tum and position. In this case, the Wigner function [8j|9] 
can be employed as a quantum generalization of a classi- 
cal probability distribution. This function is a map be- 
tween phase-space functions and density matrices. Even 
if the Wigner function can take on negative values, by 
integrating out either the position or the momentum de- 
grees of freedom, one obtains a bona fide probability dis- 
tribution for the conjugated variables. From this point of 
view, the Wigner function corresponding to a quantum 
state can be regarded as a quasi-probability distribution 
and interpreted as a joint probability density in the phase 
space [TO] . 

In this paper we apply quantum symplectic tomogra- 
phy to the investigation of open quantum systems [111 112] 



which, due to the coupling to an environment (bath), 
undergo a non-unitary dynamical evolution. A complete 
microscopic description of system-plus-bath dynamics is 
a complex many-body problem. Hence, as in general one 
aims at describing the dynamics of the system, only basic 
information about the bath is retained, according to the 
so-called open system approach. The state of the system 
is then expressed by means of a reduced density matrix, 
obtained from the total density matrix by tracing out the 
environmental degrees of freedom. The system dynamics 
is then governed by the so-called quantum master equa- 
tion. The master equation approach can be seen as the 
generalization of the Schrodingcr equation to the possi- 
bly incoherent evolution of a density matrix. In this case, 
the generator of the time evolution is the Liouville dis- 
sipative operator. The integration of a time-dependent 
Liouvillian being a highly involved task, e. g. see [THIRTI] . 
it is highly preferrable to deal with a time-independent 
Liouvillian, i.e. to assume a Markovian dynamics. Sev- 
eral approximations allow a Markovian description, such 
as the weak coupling limit , the singular limit and the low 
density limit [TT1 IT21 IT5] . 

Nevertheless, a proper derivation of the master equa- 
tion still requires complete information about the bath. 
The lack of this knowledge leads to the derivation of 
phenomenological master equations with unknown coef- 
ficients. Indeed, recent investigations [TBJ[T7J[T5] provide 
a more accurate approximation than the weak coupling 
limit, due to a more refined coarse grained dynamics. 
Even in this case, the obtained master equation has un- 
known coefficients, as it depends phcnomcnologically on 
the system investigated. 

In this paper, we will focus on a class of Markovian 
master equations with unknown coefficients modeling a 
one-dimensional damped harmonic oscillator. In particu- 
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lar, we choose Lindblad operators [H5] linear in both mo- 
mentum and position degrees of freedom, such that the 
dynamical evolution of the system preserves the Gaussian 
form of the states. Our goal is to show how, by means of 
a tomographic approach, it is possible to measure indi- 
rectly the unknown coefficients by using Gaussian wave 
packets as a probe. 

This paper is organized as follows. In section [n] we 
introduce the class of master equations we want to in- 
vestigate. In section |III| we derive the expressions for 
the coefficients of the master equation as a function of 
the first and second evolved momenta (cumulants) of a 
Gaussian state. In section IV we introduce the Wigner 



function and the Radon transform for an arbitrary Gaus- 
sian wave packet at a generic time t . We show that in 
order to measure the cumulants of a Gaussian state, and 
then indirectly the unknown parameters of the master 
equation, we need only a finite number (eight or ten) 
of time-independent tomograms. In section [V] we sum- 
marize and discuss our results and outline some feasible 
applications. Finally, in appendix A, we propose an alter- 
native procedure to obtain the cumulants of a Gaussian 
state by means of time-dependent tomograms. This ap- 
proach however appears to be less convenient for practical 
implementations. 



II. DESCRIPTION OF THE SYSTEM 

We want to investigate a class of master equations de- 
scribing a Gaussian-shape-preserving (GSP) evolution of 
a quantum state. In the Markovian approximation, the 
non-unitary time evolution of a quantum system is de- 
scribed by the following general master equation [T9]l2"0]: 



dp(t) 
dt 



L{p{t)) 



Vfp(t)V , (i) 



wherep(i) is the reduced density operator of the system. 
Eq. (Il| is exactly solvable if the Lindblad operators Vj 
and the system Hamiltonian H are, respectively, at most 
first and second degree polynomials in position (q) and 
momentum (p) coordinates [5TJ [55] . 

For systems like a harmonic oscillator or a field mode 
in an environment of harmonic oscillators (i.e. collec- 
tive modes or a squeezed bath), H can be chosen of the 
general quadratic form 



5 , 

H = H + -(qp + pq) 



Ha 



2iif + ^ (2) 



where 5 is the strength of the bilinear term in q and p, 
m is oscillator mass, and lu its frequency. The operators 



Vj, which model the environment, are linear polynomials 
in q and p: 



Vj = ajp + bjq, j = 1, 2, 



(3) 



with cij and bj complex numbers. The sum goes from 1 to 
2 as there exist only two c-linear independent operators 
V±, V2, in the linear space of first degree polynomials in p 
and q. We can safely omit generic constant contributions 
in Vj as they do not influence the dyamics of the system. 

Given this choice of operators, the Markovian master 
equation ([l]) can be rewritten as: 



dm 

dt 



Ho,p(t) 



+ 



1 

h 

i(X - 5) 



i(X + S) 



2ti 



D 



2h 

\p,p{t)q + qp{t)] 
D, 



[q,p(t)p + pp(t)] 



^([q,[p, p (t)]] + [p,[q,p(t)}}) 



(4) 



where A = 
stant and 



D 



j=i 2 a jbj is the unknown friction con- 



<i<i 
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3=1,2 



j=l,2 



D, 
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Re a *j b 3 

j=l,2 



(5) 



are the unknown diffusion coefficients, satisfying the fol- 
lowing constraints which ensure the complete positivity 
of the time evolution EH [22] : 



i)D qq > 0, ii)D pp > 



Hi^JDqqJDpp 



D 



IP 



> A 2 fi 2 /4. 
(6) 

Markovian GSP Master equations of the form Eq. Q 
are used in quantum optics and nuclear physics |23[ 124) 
I25j . and in the limit of vanishing w can be employed for a 
phenomenological description of quantum Brownian mo- 
tion [2l)l [2T1 155] , Also, in the case of a high-temperature 
Ohmic environment the time-dependent master equation 
derived in |13l 114] can be recast in this time-independent 
shape. It must be noted however that in the high- 
temperature limit the third constrain in |6| seems to 



Nevertheless. 



even if D qq = 



0, Dgp = 



be violated. 

and A 7^ 0, D pp diverges only linearly with tempera- 
ture. Therefore, we can recover the complete positivity 
by means of a suitable renormalization. This renormal- 
ization consists in adding a suitable subleading term D qq 
(e.g. D qq oc T" 1 ). Otherwise, we can consider an high 
frequency cut-off for the environment |13l 114] . In this way 
the master equation is not Markovian anymore. Anyway, 
since it involves only regular functions, it should give a 
completely positive dynamics (as the microscopic unitary 
group does). 
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III. GAUSSIAN STATES EVOLUTION 

In this section we investigate the evolution of an ini- 
tial Gaussian state according to Eq. Q. In particular 
we derive invertible expressions for the cumulants of the 
state at a time t in terms of the parameters of the master 
equation. Due to the Gaussian shape preservation, the 
evolved state at time t is completely determined by its 
first and second order momenta: 



(4)t 


= TrG3(t)g), 




(Ph 


= Tr(p(t)p), 




Aq 2 


= Mp(t)f)- 


- %l 




= Tr(p(t)P 2 )- 


- (p)l 


(q,p)t 


= Tr(p(*r 


+ pq\ 
2 ) 



-(q)t(P)t. 



(7) 



Due to the linearity of the Vj's in phase-space, the time- 
evolution of the first and second order cumulants can 
be decoupled. We then obtain the following two sets of 
solvable equations [2U [22] : 



(q) t = -(X-6)(q) t + ±(p) t 



±{p)t =-rmo 2 (q) t -(X + S)(p) t 



(8) 



fM = 



-2(A 
-2(A + 5)Ap 2 



6)Aq 2 + -a(q,q) t 
m 



2D, 



-muj^Aqf H Apt 

m 



2muj 2 a(q,p) t + 
2 2X<j(q 7 p) t 



•i 'i 



2D 



pp 
2D 



QP 
(9) 

The above equations allow to obtain the time-dependent 
momenta as a function of the Master equation coefficients 
A, D qq , D pp , D qp . We now show how to invert these rela- 
tions in order to express the parameters A, D qq , D pp , D qp 
as a function of the evolved cumulants at an arbitrary 
time. The solution of Eqs. ^ is given by [211122] 



(q)t = 



(P)t 



-A/ 



(q)o ^cosh rjt -\ — sinh rjt 

-(P>0— sinh rjt 
mrj 

mu> 2 . 

- Wo smhr/t 

V 

~{ft)o (cosh rjt sinh rjt 

V v 



(10) 



from the second set of equations (|9| , whose solutions can 
be expressed in a compact form as 



X(t) = (Te Kt T)X(0) + TK~\e Kt - l)TD, (11) 



where 



X(t) 



D 



T = 



K = 



/muAq 2 ^ 
muj 

— [ 8 - rj (5 + 77 



-2(X-rf) 
-2(A + ry) 




(2mioD qq ■ 
2D 



\ 2D qp 






-2A 



(12) 



^From the invertibility of matrices T (T 2 = 1) and K = 
K^ 1 (e Kt — l) (invertible for bounded K also if some of 
its eigenvalues are 0), we can derive the expression of 
D qq , D pp and D qp using Eq. (fill: 



D = TK~ 1 T (X{t) — (Te T)I(0)) , 
K = K- 1 (e Kt ~l) 



V 



2(X- V ) 








_ e -2(A+»,) 

2(A + r/) 




-2A 



(13) 



(14) 



2A / 



We emphasize that the time t at which we are consider- 
ing the cumulants is completely arbitrary. For instance, 
the expression of the coefficients D qq , D pp , D qp in terms 
of the asymptotic second cumulants and the parameter 
A reads: 



D„ 



(A — S)Aq l 



1 



&(q,p)c 



Dpp = {X + 5)Apl +mu} 2 a(q,p) c 



Dqp 



IV. 



mu) Aq 



m 



-Ap 2 O0 +2Xa(q,p) c 



(15) 



CUMULANTS RECONSTRUCTION 
THROUGH TOMOGRAPHY 



where rj 2 = S 2 



If rj 2 < we can set rj = iQ, and the 



previous equations hold again with trigonometric instead 
of hyperbolic functions. The coefficient A can then be 
obained by inverting Eqs. (10). 



The elements of the diffusion matrix can be retrieved 



In this section we introduce a procedure based on sym- 
plectic tomography in order to measure the first and sec- 
ond cumulants of a Gaussian wave packet at an arbi- 
trary time t. This will allow us to indirectly measure 
the parameters A, D qq , D pp , D qp , them being functions of 
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the evolved cumulants at an arbitrary time (see previous 
section). The tomographic approach is very useful when 
dealing with a phcnomcnological master equation of the 
form of Eq. Q , as the dependence of the coefficients of 
the master equation from the physical parameters is in 
principle unknown. 



A. Symplectic tomography 

Given a quantum state p{t) its Wigner function reads: 
W(q,p,t) 



1 



+°° / i2 py N 

dyexp ( — ^- ) p(q-y,q + y,t). 



(16) 



If the system dynamics is described by the master equa- 
tion Q, and the initial state is Gaussian, the Wigner 
function preserves the Gaussian form of the state. In- 
deed, it can be expressed as a function of its first and 
second order momenta: 

W(q,p,t) 



x exp 



2ir ^Aq 2 Ap 2 -a(q,p) 2 

Aq 2 t {p~{p) t ) 2 + Ap 2 {q-{q) t ) 2 
2[Ag?Ap? - a(q,p) 2 } 

2<j(q,p) t (q-(q) t )(p-(p) t ) 
2[Aq 2 Ap 2 ~a(q,p) 2 } 

Let us now consider the line in phase-space 
X — pq — vp = 0. 



(17) 



(18) 



The tomographic map of a generic state along this line, 
i.e. its Radon transform, is given by: 

w(X,p,v) = {8{X-pq-vp)) 

W(q,p, t)S {X — pq — vp) dqdp. 



JR 2 



(19) 



^From equation (|17|) it follows that for a Gaussian wave 
acket the Re 

w(X, p, v) 



packet the Radon transform can be explicitly written as: 

1 



/2ny/Aqfp 



2, ,2 



Ap\v 2 



2v(q,p)tfJ-v 



exp 



(X - p(q) t - v{p) t f 



(20) 



2[Aq 2 p 2 + Ap$v 2 + 2a{q,p) t pu] 

with the following constraint on the second cumulants: 

Aq 2 p 2 + Ap 2 v 2 + 2<r(q, p) t pv>0. (21) 

This constrain is obeyed for each value of the parame- 
ters p and v iff Aq 2 Ap 2 — <r(q,p) 2 > 0. This inequality 
is always satisfied as a consequence of the Robertson- 
Schrodinger relation. 



Eq. (|19| also implies a homogeneity condition on the 
tomographic map: |c| vj{cX, cp, cv) = zu(X, p, v). This 
condition can be used in the choice of parameters p, v. In 
fact, if one uses polar coordinates (r, 6), i.e. p — rcosd, 
v = rsin(9, the homogeneity condition can be used to 
eliminate the parameter r. From Eq. (l8| it emerges 
that the coordinates of the phase space need to be prop- 
erly rescaled in order to have the same dimensions. For 

instance, we can set q — > y/^jf-q and p — > \J nmu P' ^ n 
particular if to = 0, i.e. for a free particle interacting with 
the environment, we can choose the same rescaling with 
a fictitious frequency defined by hu) — Ap\j2ra, imposing 

-p. In general, every rescaling 



i $i« and p 



assigning the same dimensions to q and p is suitable for 
our purpose. 



B. Prom tomograms to cumulants 

Let us now consider the tomograms corresponding to 
two different directions in phase space, i.e. to two dif- 
ferent couples of parameters (p,v), e. g. X = q and 
X = p. These lines in phase space are associated respec- 
tively to the position and momentum probability distri- 
bution functions: 



zu(X,l,0) 
tn(X,0,l) 



1 



Aq t V2ir 
1 



A Pt V2TT 



exp 



exp 



2Aq 2 

[X-{p)tf 
2Ap 2 



(22) 
(23) 



^From Eqs. ( 22 1- ( 23 1 we see that the tomographic map 
depends only on a single parameter X. This reduces 
the dimensionality of the problem with respect to the 
Wigner function, that is a function of both p and q. 
The lines individuated by the choices (p.v) = (1.0) and 
(/i, v) — (0, 1) correspond to tomograms depending on 
the time average and variance respectively of position 
and momentum. In order to determine the latter quan- 
tities we have to invert Eq. ( 22 1 and ( 23 1 for different 
values of X, i.e. for a given number of points to measure 
along a tomogram. Thus, our first goal is to determine 
the number of tomograms required to measure the cumu- 
lants of our Gaussian state. 

To answer this question, we first focus on the direction 
p = 1, v = 0. In Fig. [I] we plot the Wigner function 
of our system at a generic time t and some straight lines 
along the considered direction. In Fig. [2] we plot the GSP 
tomogram defined by Eq. (22 I. Inverting Eq. (22 I, we 
obtain: 



(X-(q) t Y 



2Aq 2 t In 



vj{X, 1, 0) Aq t V2tt 



(24) 



Using the value of the tomogram tn(0,l,0) we can get 



5 



(q)t as a function of Aq t : 



(q) t = ±Aftj21n 



tn(0, 1, 0)A gtV / 2^ 



(25) 



If we know the sign of (q)t then we need only the value 
of the tomogram tu(0, 1, 0) to get (q) t , otherwise we need 
another point. Using Eq. (25), Eq. (22 1 becomes an equa- 
tion for Aq t only, and it can be rewritten as 



2Aq^ In 



1 



w{X,l,0)Aq t V*r 



XT Ag M /21n 



ro(0,l,0)A ?tV / 27f 



(26) 



This equation is trascendental, therefore we will solve it 
numerically. We can graphically note in Fig. [3] that for 
each X and corresponding w[X, 1,0) there may be two 
values of Aq t satisfying the previous equation. In order to 
identify one of the two solutions, it is enough to consider 
two points, {(Xi,w(Xi,l,0))} and {(X 2> vj{X 2 , 1, 0))}, 
and to choose the common solution for the variance. This 
is made clear by Fi g. [3] where the ratio between right 
and left side of Eq. ( 26 1 for two different values of X is 



plotted. The common solution (i.e. when both ratios are 
equal to 1) is labeled Aq t . 



W(q,p,t) 




f 2Ap, 



p-(P)t 



q~{q)t 



FIG. 1: Wigner function, W(q,p,t), of Eq. jl7[ l and some straight 
lines along the direction fi = 1 and v = on the plane qp. 
Ag 2 Ap t 2 - a(q,p)l = 0.64Ag 2 Ap 2 . W M = , .\ ; =■ 

27r V A lt ~a(q,p)i 



As a consequence, whether we know or not the sign of 
the average (q)t, we need three or four points to deter- 
(q)t and Aq t in Eq. 



mine [q) t and l\q t m iUq. (22 1. Analogously, we need 
other three or four points for(p) t and Ap t in Eq. (23 1. 
Let us now compute the covariance a(q,p)t- To this 



w{X,l,0) 




-Lq, kq 



FIG. 2: Tomogram, m(X, 1,0), of Eq. |22| for the direction = 1 
and v = 0. 




FIG. 3: The ratio between right and left side of Eq. l ]2(j[ l for two 
different values of X with = 1 and v = is plotted as a function 
of Aqt. Values used: (q)t =3 and Xi = 4.5 (continuous line) and 
X2 = 2.5 (dashed line). The values of 1,0) are computed 

using Aq t = 1 to simulate what one would get experimentally. 
Continuous line shows how Eq. ( |26| with X\ = 4.5 and ro(4.5, 1, 0) 
can be satisfied (ratio=l) by two values of Aqt- The comparison 
with a second case with X2 = 2.5 allows to determine which values 
of Aqt is the right one. 



purpose, we consider the tomogram: 



2V2/ W A 9? + Ap 2 + 2a(g,p) t 



exp 



Ag t 2 + Ap 2 + 2(7(9, p) t 



(27) 



This is a Gaussian whose average value is already de- 
termined. Indeed, according to the previous steps, we 
need two more points of this tomograms to determine 
the spread (Aqf + Apj)/2 + a(q,p) t , from which we can 
retrieve a(q,p) t - 

Hence, we have shown that by means of eight or at 
most ten points belonging to three tomograms, the first 
and second order momenta of a Gaussian state can be 
measured at an arbitrary time t. One can then use these 
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measured cumulants in order to infer the master equa- 
tion parameters describing the system under investiga- 
tion. We note also that we can reasonably infer that 
the number of tomograms needed to reconstruct the sys- 
tem density operator is minimized by employing Gaus- 
sian wave packets as a probe. Indeed these states have 
minimum uncertainty, and are the only states having pos- 
itive Wigner function [25] , 



can re-write it as a combination of four Gaussian func- 
tions. Therefore, due to the linearity of the master equa- 
tion, it can be possible to derive analytically the state 
evolution and to analyze its loss of coherence by means 
of the procedure we propose. 
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V. CONCLUSIONS 

In this paper we have proposed an approach to the 
study of open quantum systems based on quantum sym- 
plectic tomography. 

In many contexts the reduced dynamics of a system 
coupled with its environment is modeled by phenomeno- 
logical master equations with some general features, but 
with unknown parameters. Hence, it would be highly ap- 
pealing to find a way to assign some values to these pa- 
rameters. We have tackled this problem for a wide class 
of Markovian master equations, which are the Gaussian- 
shape-preserving ones. We have proved that it is possi- 
ble to retrieve their unknown parameters by performing 
a limited number (ten at maximum) of time-independent 
measurements using Gaussian wave packets as a probe. 

This result leads to some interesting applications. 
Once retrieved the unknown master equation coefficients, 
it is possible to compute the dynamical evolution of any 
physical quantity whose analytical expression is known. 
The indirect-measurement scheme we propose could be 
then employed to make predictions on system loss of co- 
herence due to the external environment. In order to per- 
form this kind of analysis one can consider some quanti- 
ties such as the spread and the coherence length in both 
position and momentum |30) . provided their analytical 
expressions are available for an arbitrary time t (e.g. see 
Ref.|28j). Working in the coherent state representation, 
the evolution of the system of interest from an arbitrary 
initial state can be in principle predicted. Therefore, it 
is possible to perform the proposed indirect analysis of 
the decoherence processes. For example, if we consider 
an initial Schrodinger-cat state, highly interesting due to 
its potentially long-range coherence properties and its ex- 
treme sensitivity to environmental decoherence |31j . we 
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APPENDIX A: ALTERNATIVE PROCEDURE 

Here we propose an alternative time-dependent pro- 
cedure to compute the second cumulants of a Gaussian 
state, by means of tomograms, given the knowledge of 
the first cumulants time evolution. To this purpose we 
need to consider the following three tomograms: 



ZD I 
— 2 
G7 3 



U7((p) t) Q,l) 

^((<?>t,l,0) 



1 



1 



/2nAq t 

( {ph + Wh ± 

V \/2 ' \/2' \/2y 
1 

2Hy/Aq?/2 + + <j(q,p) t ' 



(Al) 



Inverting the previous equations one can infer Aq t , Ap t 
and a(q,p)t from the knowledge of w\, wi and 7373. How- 
ever, this procedure presents two drawbacks. In fact, the 
evolved averaged values (q) t and (p) t are required and we 
need tomograms evaluated on time-dependent variables. 
These problems do not arise in the time-independent pro- 
cedure, based only on tomograms for which no a priori 
knowledge on the Gaussian state is required. Neverthe- 
less, in this alternative time-dependent scheme only three 
tomograms are required. 
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